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Abstract 

Isosinglet fermions naturally arise in a variety of extensions of the Standard Model, in 
particular in models with extra dimensions. In this paper, we study the effect of the ad- 
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ton flavour violation in neutral current and lepton universality constraints are extensively 
discussed. We show that new physics effects in v e — v„ CP asymmetries are significantly 
enhanced due to leptonic maximal mixings but still too small to give a signature at future 
neutrino factories. A signal for CP asymmetries in z/„ — v T channel due to new physics 
could be observed at 1 — 3 a if lepton flavour violating r decays are seen in a very close 
future in U-factories like BELLE experiment. 
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I. INTRODUCTION 



Two of the most exciting recent developments in particle physics are, on one 
hand, the experimental evidence for neutrino oscillations [1, 2] pointing towards non- 
vanishing neutrino masses and mixings and, on the other hand, the suggestion of 
models with more than three spatial dimensions [3] which have the attractive feature 
of addressing, within a novel framework, long standing puzzles like, for example, the 
gauge hierarchy problem and the elementary fermion mass spectrum. 

In this paper, we analyse some of the phenomenological consequences of isosinglet 
charged leptons which naturally arise in some of the models with extra-dimensions, 
discussing in particular the interplay between the study of lepton flavour violating 
rare processes and new physics effects observable in CP asymmetries at neutrino 
factories [4]. 

In models with extra-dimensions, the use of new techniques is accompanied by 
a very interesting phenomenology, with a large number of new particles populating 
the spectrum at energies typically of the order of the compactification scale. Extra- 
dimensional models represent then a very interesting framework for physics beyond 
the SM, with the Kaluza-Klein (KK) excitations of the bulk fields being extra gauge 
bosons, fermions or scalars. 

Arkani-Hamed and Schmaltz [5] have proposed in the context of these theories 
the split fermion idea, which stands for the use of the localization properties of bulk 
fermion zero modes to suppress coefficients without the need of symmetries protecting 
the corresponding operators. One particularly interesting application of this idea is a 
natural realisation of the observed pattern of fermion masses and mixing angles [6, 7]. 
Recently, it was shown in [7] that, under very general circumstances, it is possible to 
decouple the mass scale of the fermion KK excitations from the compactification scale 
by means of multi-localizing the fermion zero modes. In this way it was possible to 
naturally realise the observed spectrum of quark masses and mixing angles in a model 
with a multi-brane scalar background, a compactification scale M c ~ 100 TeV 1 and 
still have a light KK excitation of the right-handed (RH) top quark with observable 
phenomenological consequences at present or future colliders. 

1 This value is consistent with the bounds on flavour changing neutral currents in these models [8] . 
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In the class of models described in [7] , the RH component of the heaviest fermions 
are naturally multi-localized. Thus they typically have light (vector-like) KK exci- 
tations mixing mainly with them. One extra ingredient appearing in these models 
is that the SM fermion gauge couplings are modified due to mixing with their KK 
excitations, with corrections which are proportional to the masses of the SM fermions. 
We will apply the above scenario to the leptonic sector [6, 7], analysing in detail the 
phenomenological consequences of the isosinglet charged leptons which naturally arise 
in this class of models. The most natural situation is to have a multi-localized RH 
tau lepton with a light isosinglet KK excitation mixing mainly with it. We will see 
in the sequel that, for phenomenological reasons, it is crucial that the new vector-like 
charged lepton mixes mainly with the tau. 

The analysis of new physics effects in the leptonic sector, is specially relevant in 
view of the various experimental projects designed to uncover the neutrino world, 
in particular the planned neutrino factories [4] which will measure neutrino masses 
and mixings with great precision and eventually CP violation in the leptonic sector. 
Furthermore, neutrino factories have the potential to uncover new physics, which 
can arise in a variety of extensions of the Standard Model (SM), such as models with 
vector-like neutrinos [9]. In the class of models with vector-like charged leptons which 
we consider in this paper, we will find explicit realizations of the model independent 
results previously found [10, 11]. 

The paper is organised as follows. In the next section, we present a general de- 
scription of models with extra isosinglet charged leptons, with special emphasis on 
lepton-flavour-violating processes and introduce the structure motivated by models 
with extra dimensions. The results of this section are model independent and there- 
fore relevant for any extension of the SM with isosinglet charged leptons. In section 
III, we discuss the phenomenological implications of isosinglet charged leptons on 
neutrino oscillations and in particular on CP asymmetries at neutrino factories. In 
section IV, we present our main results while our conclusions are summarised in sec- 
tion V. We end the paper with an appendix where the KK description of theories 
with extra dimensions is reviewed. 
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II. GENERAL DESCRIPTION 



In this section, first we recall some well-known features of models with vector-like 
particles, in particular the structure of the neutral and charged flavour-changing lep- 
tonic currents. Secondly, the constraints coming from rare lepton-flavour violating 
decays and lepton universality are presented. Finally, the extra-dimensional frame- 
work for isosinglet charged lepton is described. 

A. Charged lepton masses and mixings. 

Let us consider the addition of a vector-like charged lepton, singlet under SU(2)l, 
to the SM spectrum. The mass Lagrangian for the charged leptons can be written 
with complete generality in the current eigenstate basis as 



where mi is a 3 x 3 mass matrix and J describes the mixing between the first 3 families 
and the vector-like charged lepton. Ll,r are the Left-Handed (LH) and Right-Handed 
(RH) components of the isosinglet charged lepton and l' L , Ir the usual LH and RH 
charged lepton fields respectively (e,/i,r). Note that the (1 x 3) zero matrix in 
eq.(l) corresponds to an allowed choice of weak-basis and does not imply any loss of 
generality. The matrices mi and J are A J = 1/2 mass terms, therefore proportional 
to the SM Higgs vacuum expectation value while M is an SU(2)l ® U(l)y invariant 
mass terms, which can be arbitrarily large, since it is not protected by the low energy 
gauge symmetry. 

Before diagonalising the charged lepton mass matrix, let us have a look at the 
neutrino mass sector. We shall assume that naturally small LH neutrino Majorana 
masses are generated through the breaking of lepton number at high energy. It is well 
known that the seesaw mechanism provides one of the most attractive scenarios for 
generating small LH neutrino masses [12]. However, for our discussion, the detailed 
origin of neutrino masses is not important. The Majorana mass term for the 3 light 




+ /i.e., 
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LH neutrinos can be written as 

£„ L=2 = \v c lM„vl + h.c, (2) 

with M v a 3 x 3 symmetric matrix which can be diagonalised by a unitary transfor- 
mation U v 

M V U V = diag(m Va ,m v ^rn VT ). (3) 

We can go to the mass eigenstate basis for the LH neutrinos by making the corre- 
sponding unitary transformation on the lepton doublets 



h 




(4) 



Of course, eq.(l) is not invariant under this transformation. In the neutrino mass 
eigenstate basis, the mass matrix for the charged leptons is given now by 

C m = ( l L L L ) ( ^ ^ j ( ^ j + hx - = ChWfa). (5) 

It may seem not necessary to discuss the neutrino sector at this point. Indeed, without 
loss of generality, we could have started our discussion with eq. (1) in the physical 
basis for neutrinos. However, models where vector-like fermions naturally appear 
(e.g. models with extra-dimensions) often predict some textures for J, mi or M v . In 
such cases, it is convenient to work with eq. (5) where the constraints on the structure 
of J, mi and M v can be implemented in a straigthforward way. 

The charged lepton mass matrix M\ is diagonalized by the unitary transformations 
Ul and Ur 

U\M l U R = diag (m e ,m /J ,m T ,M D ), (6) 
with the mass eigenstates 1™ R given by 

(Il,r) = U L , R {IZ R ). (7) 
The left rotation Ul can be written as [13, 14] 

UlK Up/M 



I -,PKIM 1 / 



-.Pk/m 

where K is the unitary matrix which diagonalizes m\m\. 



B. Z-mediated flavour changing neutral current interactions. 



The leptonic neutral current gauge interaction reads, in the weak eigenstate basis, 



£z = -V^(^ - sin 2 ^JL), (9) 

COS (7,„ 



with 



= ~ Ih-flL, (10) 

jl = -(m + L^L). (ii) 

In the mass eigenstate basis, for the light charged leptons, one gets 

£z 9ht = -tt^Z&UI - 2 sin 2 9 W ) - & k ]ri? k , (12) 

Z COS u w 

with i, k — 1, 2, 3, (3ik = U^Ulm- The effect of the vector-like singlet is the appear- 
ance of the flavour changing neutral couplings (FCNC) flu- at tree level. Note that at 
first order in mf/M 2 , 

a _ J l K U J k K kj /-i q\ 

j ~ ' { 1 

and thus these couplings are naturally suppressed if J/ -C M. 



C. Higgs mediated FCNC. 



The interaction between charged leptons and the neutral Higgs boson is given, in 
the weak eigenstate basis, by 

9 n ... , , 7 rr , l -1 rrO 



c 
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[l L mil R + l L JL R + h.c] H c 



(14) 



2M W 

After diagonalisation of the mass matrix, the Lagrangian for the light charged leptons 
becomes 



p light 
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2M, 



w 



(15) 



I Li ( m ihk ~ MUl M U mk ) Q k + h.c. 
where i,k — 1, 2, 3 and m 12i 3 = m ettJi)T . The interaction with the pseudoscalar neutral 
field x is given by 

[l L mil R + L L Jl R - h.c] x- 



(16) 



2M W 

Proceeding in the same way as for the neutral Higgs scalar, the interaction with the 
pseudoscalar Higgs field x is given by 



n light 



2M 



w 



l Li (mAk - MU* LAj U mk ) Z™ - h.c. 



(17) 



D. Flavour changing charged current interactions. 



The interaction with the W is given as usual by 

C w = ^=(W-J^ + h.c), (18) 

with 

j^+=7 L yv L . (19) 

In the mass eigenstate basis, the charged current interaction reads 

Cw = ^( W X^7X* + h.c), (20) 

with a = 1 ... 4 and k = 1 . . . 3. The Vmns mixing matrix is therefore 4x3 and given 
by 




E. Limits from rare lepton flavour changing decays and lepton universality. 

Lepton-flavour violating processes are strongly constrained by the experimental 
limits on rare /x and r decays (see for instante Table I). In the present model, Z- 
mediated FCNC induce tree level corrections to processes like \i — > 3e or \i — > e in 
^Ti while \i — > ej receives corrections at one-loop level. We shall consider first the 
limits coming from tree-level rare tau and muon flavour changing decays and from 
gauge coupling lepton universality. Limits coming from loop-induced processes like 
H — > e7 will be discussed afterwards. 

TABLE I: Experimental limits on rare \i and r decays constraining FCNC. 



Br(fjL - 


+ erf) < 1.2 x 10" 


-ii 


[15] 


Br(fi - 


+ 3e) < 1.0 x 10- 


-12 


[16] 


Br(n — ► e 


in ^Ti) < 6.1 x 


10" 


[17] 


Br[r - 


-» 3e) < 7.8 x 10 


-7 


[18] 


Br(r - 


-> 3/i) < 8.7 x 10 


-7 


[18] 
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1. t and [X rare flavour changing decays and lepton universality 



(23) 



The restrictions on rare lepton decays can be translated to indirect bounds on lep- 
ton flavour violating branching ratios for Z decays [19], which are directly connected 
to the parameters of the Lagrangian. For instance the branching ratio for Z — > Ulj is 
approximately given, for i ^ j, by 

BriZ^kl,)^ 1 -^ 2 . (22) 
Using the current experimental limits, one obtains 

I A* | < 1-1 x 1(T 6 from (// -> 3e), 

|/^ e | < 4.0 x 1(T 7 from (p -»■ e in <gTi), 

\(3 Te \ < 2.3 x 1(T 3 from (r -> 3e), 

\Prr\ < 2.4 x 10" 3 from (r -> 3/i). 

Lepton universality sets limits on the diagonal elements f3„. Indirect bounds from 
unitarity loss in charged currents have been computed for the case of vector-like 
neutrinos [9] but give much less stringent constraints than direct violations of univer- 
sality appearing in the model with extra isosinglet charged leptons. Using the one 
sigma deviations for the effective charged lepton couplings, gy 1 ^ 7 , to the Z [20] as an 
estimation of the allowed new physics contribution we get 

\Pee\ < 0.0007, 

\Pw\ ^ °- 0011 ' ( 24 ) 

\P TT \ < 0.0013. 

The last two limits are restrictive enough to impose a bound on the f3^ T coefficient 
stronger than the one obtained from rare tau decays 



|/V| = y|A^|/3 TT |<1.2xl0- 3 . (25) 

These limits can be translated to constraints on the elements of the charged lep- 
ton mass matrix making further assumptions. As examples, we shall consider two 
theoretically motivated ansatze for Jj. 

Case A 
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Given the fact that experimental constraints are much more stringent for the first 
families than for the heavier ones and that both the entries in the mass matrix J and 
the standard lepton masses arise from Yukawa couplings, a reasonable ansatz for the 
values of the J\ is [21] 



^ ~ ^i. (26) 

Ji m 

Using eq.(26) and assuming that the mixing induced by K is negligible, it is possible 
to extract limits on Ji/M from the experimental bounds given in eqs. (23,24), the 
strongest bounds being 

Si = ^- < 4.4 x 10~ 5 , 



M 

5 2 = < 0.91 x 1(T 2 , (27) 

5 3 = -4 < 0.036. 

These constraints come from /x — > e conversion in nuclei for the first two families and 
from lepton universality for the third one. 

Case B 

It is however possible to use a different ansatz, motivated by models with extra 
dimensions [7, 22] reviewed in the Appendix, which leads to a natural suppression for 
the corrections to the couplings of the first family. Let us assume that in the weak 
basis, the vectors J, and (m;) i3 are parallel, 

Ji = A(m,) l3 , (28) 

with A a flavour independent constant. In this situation we can use the following 
property to compute the value for the (3^ parameters, 

(Jttf)* = 4-J, = A4(m ( ) 33 = X mi (K R )l 3 + O Q) 2 , (29) 

with no summation on the index i and where m ; are the mass eigenvalues of the 
charged leptons and K R is defined by 

K ] miK R = diag(m ei m M , m T ) + O (J^pj . (30) 

Then, we can write 

Pij = A 2 ~ } M \ 3J 3 - (31) 



The ratio of FCNC processes is then not only suppressed by masses but also by mixing 
angles. For instance, 

(32) 



I&A.I rn e \Kg\ 



\M m T \Kgy 

and, apart from the mass suppression, there is an extra suppression by mixing angles 
ii\Kg\^\K*\. 



2. fi — > ey 

This decay is induced at one-loop level, the constraints coming from this process 
are thus expected to be less stringent than the ones arising from the processes studied 
in the previous section. Nevertheless, the experimental progress on \x — > expected 
in the next few years [23], makes it worth studying it in detail. The branching ratio 
for the transition fi — > has the form 

4 

Br{fi -> ej) = 3847r 3 a^- |A(m 2 )| 2 , (33) 

m fj, 

where v = (SG^)" 1 / 4 - 174 GeV is the Hi gg s 

vacuum expectation value, a is the fine 
structure constant and A(q 2 ) the form factor coming from the one- loop computation. 
At one-loop order, the transition amplitude is induced by 3 differents classes of Feyn- 
mann diagrams respectively due to W, Z and Higgs exchange. For diagrams with W 
exchange, the photon line is attached to the W's. This diagram is suppressed in the 
SM by neutrino masses whereas in our model, due to the non-unitarity of the Vmns, 
it gives an independent contribution. Let us consider the form factors arising from 
these three kinds of diagrams 2 

A = A w + A z + A H - (34) 

We evaluate these form factors in the limit w? lv jm^ v z H —>■ for the light families 
and we assume that in our model Mo 3> mw,z,H for the isosinglet charged lepton 
contribution 3 . They read, in this limit, 

2 See ref.[24] for details on the form factor computation. 

3 Indeed, experimental search on heavy charged lepton imposes > 100.8 GeV [25]. 
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a * - (I - IS) • < 36 > 



The branching ratio for \x — > &y is then given by 



B ^ - erf) ~ — l^^l 2 (^24 + T^dr + Zt£ ) 



48 Ml, 6M2 



Using the strongest limit on \UL4^,Ul 4e \ implied by eqs.(23) and the fact that M|,/m| H 
> 1, this branching ratio is bounded by 

Br(fi -> e7) < 5.5 x 1(T 16 , (39) 

far below the sensitivity of current and planned experiments. As expected, fi — > e'j 
does not introduce any significant restriction, in a model with isosinglet charged 
leptons. 



F. Extra-dimensional framework for isosinglet charged leptons 

As we have emphasized in the introduction, theories with extra dimensions pro- 
vide an exciting arena for the study of physics beyond the SM. It has been recently 
shown that models with multi-brane backgrounds [7] naturally present light isosin- 
glet fermions which correspond to the first KK excitations of the heavier SM weak 
singlets (tr in the leptonic sector). We review in the appendix the details of the 
KK description of the theory, the relevant features at this point being the low energy 
spectrum and the form of the Yukawa couplings. In the case of interest, the spectrum 
consists on the fermion zero modes, which are the SM fermions, plus one relatively 
light (as compared to the compactification scale) isosinglet charged lepton, mixing 
mainly with the tau. The rest of the spectrum have masses of the order or the com- 
pactification scale which decouples from the low energy physics. The other feature 
relevant for our discussion is the form of the Yukawa couplings. We take the scalar 
field responsible for the electroweak symmetry breaking to live in a four- dimensional 
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boundary. The effective Yukawa coupling for the n— th and m— th modes satisfy the 
following factorisation property 

,(5)e 
TlR 



\t m) = ^f l L {n) Wf e / m \0), (40) 



where is the five-dimensional Yukawa coupling and, for definiteness, we have 
taken the Higgs field to live in the y = boundary. Notice that this factorisation 
property implies the relation in eq.(28), with 

A - S^M ,41) 

As we have seen above, this property has a very relevant feature from the phe- 
nomenological point of view, which is the fact that corrections on the fermion zero 
mode gauge couplings due to the mixing with the heavy vector-like excitations scale 
with the masses of the fermion zero modes and are thus only relevant for the heavy 
fermions [22]. 

It is now clear how the hierarchical pattern of charged lepton masses is generated. 
Starting with five-dimensional (thus dimensionful) Yukawa couplings of natural or- 
der, \\j ~ ttR, we can generate effective Yukawa couplings of the appropriate size 
modifying the localization properties of the different fields and thus their overlapping 
with the Higgs boundary. Given the relatively large mass of the tau lepton it is nat- 
ural to have its RH component strongly localized near the Higgs boundary and thus 
multi-localized (see the Appendix) in the set-up we are considering 4 . 

If we ask for the RH tau lepton to be strongly multi-localized, so that it will have 
a light KK excitation, then it is natural to have all mixing angles in the charged 
leptonic sector small but the one in the (2, 3) LH sector, which can vary, depending 
on the specific values of the parameters, from moderately small to maximal. In the 
first case all bi-maximal mixing has to be generated from the neutrino sector (what 
can be easily accomplished in the model we have [26]) whereas in the second one only 
the (1, 2) large mixing comes from the neutrino sector. 



4 Notice that from the tiny neutrino masses we could expect the third family doublet not to be too 
strongly localized at the Higgs boundary, thus the mass of the tau being essentially generated by 
the order one overlapping of its RH component. 
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III. EFFECTS OF VECTOR-LIKE CHARGED LEPTONS ON v OSCILLA- 
TIONS 



A. General parametrisation 

The new interactions generated by the extra vector-like charged lepton have im- 
portant effects on neutrino oscillations. Similar effects occur in the charged current 
sector in models with extra isosinglet neutrinos [9, 27] while neutral current processes 
could allow to discriminate between them. In this section we compute the effects of 
the vector-like charged leptons on neutrino oscillations. 

As we have seen in the previous sections, in the case of one vector-like charged 
lepton, the leptonic mixing in charged currents is described by a 4 x 3 matrix Vmns- 
The mixing among the light families is given by a non-unitary 3x3 submatrix of Vmns 
which, using that deviations from unitarity are necessarily small, we will separate in 
two terms 

(V MN s)ik = {U MNS )ik + (SU) ik , i, k = 1, 2, 3, (42) 

where Umns is a 3 x 3 unitary matrix which accounts, at leading order, for the mixing 
among the light families [28] and 5U describes the small deviations from unitarity. In 
the following we consider the mass eigenstate basis for the charged leptons and the 
flavour basis for neutrinos, defined as 

v{ = U M ns<- (43) 
The charged current interaction reads in this basis, for the light leptons, 



light 9 



w-r Ll (s ik + y»L + h.c] ee c™ + c% p , m 



' w "V2L' 

where Cyy 1 , the term proportional to the identity, is equivalent to the SM charge 
current interaction and C^f describes the New Physics effects due to the isosinglet 
charged lepton with coefficient Af k p = SUij(U MNS )jk. 

In order to compute the effect of New Physics, we need to know Umns an d SU. For 
that, we have to diagonalise the charged lepton mass matrix given in eq.(5). To define 
the SM families, we shall proceed in 2 steps: first, Mi is diagonalised in blocks so that 
the heavy charged lepton mass eigenstate is defined. In the second step we make a 
unitary transformation involving only the light charged leptons, which completes the 
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diagonalization. This stepwise diagonalisation corresponds to making the following 
sequence of rotations, for the LH sector, 



Ul = -^34-^24-^14-^12-^13-^23' 



(45) 



R 



13 



Ri 



RiMj = °)- 



(46) 



where 

( ci 3 s l3 e- iSl!i \ 

10 

-s 13 e i5 ^ cia 
\0 1 / 

An explicit form for Ul can be found in ref . [29] . It is easy to check that making the 
rotation in the order defined in eq.(45), we have 



v\m x m\u l 



R^Rx^RyiR^R^^R-mMi M t -R34 i?24-Rl4 R±2 R-13 -^23 



R23R13R12 



mf 







Ml 



- fl 12- rl 13- n '23- 



(47) 



Thus -R23-R13-R12 is the mixing matrix connecting the 3 light families. It is important 
to recall that since M D 3> J 3> m T ^ e , the mixing angles 9 i4: (i = 1,2,3) are well 
approximated by 

e i4 e iS * ~ (48) 



M 

With this prescription, the Umns matrix as defined by eq.(42) is given by 



MNSJik 



(R23R13R 



12)ik, 



(49) 



and the New Physics contribution by 



-4 



NP 



MNSljk 



UmNsRuRIa -^34 U MNS 



Using the above parametrisation, Umns can be expressed as 



/ 



U 



MNS 



C12C13 



C13S12 



-C23S12 - Ci2Si 3 S2 3 e J<513 C12C23 - Si2Si 3 S23e J<513 



si 3 e 
C13S23 



<5l3 \ 



y -ci2C23Si 3 e tSl3 + S 12 S 2 3 -C12S23 - C 2 3Si2Si 3 e 4<513 c 13 c 23 J 



(50) 



P, (51) 



where P = diag(l, e 4/3l 5 e*^ 2 ). The phases flip are Majorana-type phases which will not 
play any role in our discussion. The New Physics contribution can be also explicitly 
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written as 



A NP = U 



MNS 



I cu -si 4 s 24 e- i ( <5l4 - <524 ) -s M s u c 2 Ae- i514 ^ 

c 24 -s 2 iS M e- i5 ^ 
y c 34 j 



with / denoting the 3x3 identity matrix. 

Note that the full 4x3 Vmns matrix contains, apart from 5 13 , two extra Dirac- 
type CP violating phases. These phases appear always in combination with the small 
mixing angles 9^. The appearance of these two extra Dirac-type phases has to do 
with the fact that the 3x1 column matrix J is complex, with one of the phases 
eliminated through the rephasing of the isosinglet charged lepton field [13]. 

Once we have separated the New Physics contributions from the SM Lagrangian, 
we can use the formalism described in [10, 11] to compute the effect of this new 
interaction on neutrino asymmetries. 

In neutrino factories, neutrinos are usually produced through muon decay. There- 
fore, assuming that we have a fi + beam, typically the production process of a neutrino 
state, u a , in conjunction with a e + is given by 

fi + — >■ e + u a V fl . (53) 

The production of a neutrino state u e in conjunction with a e + is also possible through 
the process 

/i + -> e + u e V a . (54) 

The detector process, for a wrong sign event, is v p d — > fi~u. Between production and 
detection, the oscillations of v a from (53) or of v e from (54) into a v p can take place. 
The important point to note is that only the process given in (53) is able to interfere 
with the SM process, fi + — > e + u e V p , due to the fact that they have the same initial 
and final states e + and V^. So in the sequel, we shall neglect the contribution 
coming from process (54) as it does not interfere with SM amplitude, i.e. its effects 
are higher order in New Physics compared to process (53). 

This approach can be generalised to any flavour and to any production and detec- 
tion processes [11]. Consider the Lagrangian responsible for the production of a up 
in conjunction with a Z+, 

C s = 2V2G F (5 aP + e s a p) (jI-fPL^) (u^P L l a ) , (55) 
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where Pl is defined as usual by (1 — 7s)/2. According our notation, 

4/3 = (A NP *U- (56) 

For the detection of a z/g signalled in the detector by a Z~, the similar 4 Fermi inter- 
action enters 

C d = 2V2G F (5 a p + e^) (u^P L d) (J a ^P L ^) • (57) 

with 

4 = 4* =► e " = ( 58 ) 
As done in ref.[10], we define v s e as the neutrino state that is produced in the 
source in conjunction with a e + , the production process being /i + — » e + t / Q ,z/ A1 , and ^ 
the neutrino state that is signalled by /i~ production in the detector, the detector 
process being v a d — > ^t~u. The detector process for a wrong sign event is interpreted 
as due to the oscillation of v e into z/ M . 



\<) = E K 1 + OOW)* + <,(tW);,, + 4r(^5)y K™) 
i 

= [V + ?)UmnsUK) , ( 59 ) 

i 

= [{I + ^UmnsUvD- (60) 

We recall that the complex conjugation comes from the fact that once Umns is defined 
in the Lagrangian as given by eq.(44), one has [30] 

W) = (u MNS y ei K) ■ (ei) 

The observation of new physics in the detection requires not only knowledge on 
neutrino beams but also on quarks properties within nuclear matter [11]. In par- 
ticular, present uncertainties on parton distribution and hadronisation for the range 
of neutrino energy will make very difficult to distinguish the effects of new lepton 
flavour violating physics from other sources of uncertainties or new physics. So, in 
the following part of the paper, we shall only compute the effects of lepton flavour 
changing interactions coming from the production processes which is purely leptonic 
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and neglect their effects in detection processes. Of course, our analysis can be trivially 
extended to include flavour changing effects in detection processes. 

The CP asymmetries in vacuum and in matter can be computed using the standard 
procedure. 

B. Oscillations in vacuum 

It is well known that in the framework of the SM, there is no CP violation in 
neutrino oscillations in the limit s\ 3 — > 0. Of course, in the presence of physics 
beyond the SM, this is not true anymore. Therefore, it is interesting to study what 
happens in our model when S13 — > 0, answering the question whether even in this 
limit, one can expect to observe CP violation in neutrino oscillations in models with 
extra isosinglet charged leptons. 

It should be recalled that FCNC impose strict constraints on the (Ul)^ = Si 
elements of the Ul lepton mixing matrix. Using the unitarity of Ul and eq.(48), we 
can relate the angles 9a to the Si coefficients 

But** * - (UlK) rj S] ~ - (ul NS ) i3 4 (62) 

The stringent bounds on new contributions to the first family couplings require S\ <C 
5*2,3 • We therefore consider the following texture for S 

S=(0,S 2 ,S 3 ). (63) 

In such a case, 9 a read 

0^*" = Sl2 (c 23 S* 2 - s 23 S* 3 ) + c 12 s 13 e~ iS ^ (s 23 S* 2 + c 23 S* 3 ) , (64) 
OyAfi^ = -ci2 (c 23 S* 2 - s 23 S* 3 ) + s 12 s 13 e- lS ^ (s 23 S* 2 + c 23 S* 3 ) , (65) 

6» t4 = - Cl3 (523^2 + °23S 3 ) • (66) 

It is interesting to note from the above equations that as expected when s 12 — > 0, 
9 e 4 S% = (Ul)^ which is strongly constrained by lepton flavour violating muon 
decays. But if s\ 2 ~ c\ 2 ~ 0(1), 9 e ^ is dominated by S 2 and S 3 contributions which 
are much less constrained than S\ (see eqs.(27)). 
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As discussed before, the model with one extra isosinglet charged lepton has two 
additional phases. For simplicity, we shall choose the phases of S 2 and £3 such that 
9 T 4 is real (£34 = 0). From these equations, the following limits can be considered 



(A) s 13 = 



9 e4 = -tan 0i2 x 6> m4 , (67) 
8u = $24- (68) 



(B) 5i3 = and s 12 = C\ 2 = 1 / %/2- We can make a perturbative expansion in s 13 

Q 

|5i4-5 24 | - 2^-s 13 sm(p + 0(s 2 13 ), 



OeiM - ±7^= =F s i3^ cos0 + 0(s? 3 )J 
sin5 14j2 4 ~ sin0 ( 1 ± Si 3 — — cos0 + 0(s^ 3 ) ] , 



(69) 



a 

with = arg(c23S , 2 * — 5235*3) and a = \c 23 S^ — s 23 Sl\. The upper and lower 
signs in the equation for 6 ei ^ correspond respectively to 6 ei and 6^ 4 . 

Using the relation between e s and A NP and expanding A NP in terms of 6^4, it is 
easy to get the expression for e s (from now on we denote Umns simply as U) 

4 = -UjxU^f - u j2 u: 2 °f - u j3 u* 3 9 -f 

-u^u^e^e 1 ^- 5 ^ - u^u^e^ - u j3 u* 2 e^e r4 e iS *\ (70) 

Following the formalism described in [10, 11], the dominant contribution to CP 
asymmetries in the different channels can be easily computed. We can use directly 
their results just replacing ef - by the expressions given in eq.(70) and efj ~ 0. 

Before studying the different channels of oscillations, let us define some notation 
that is used later, 

Am?- = m\ - rrij, (71) 



2 /yzn 

AmlL 



\j = Am?-/(2£), (72) 



AE 

Am?- L GeV 
1.27— x — x — — . 73 

eV 2 km E K J 
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To get an idea of the order of magnitude of x^-, let us evaluate them using the data on 
neutrino atmospheric [1, 2] (Amf 3 = 3 x eV 2 , tan#23 = 1) , and for Lqs — 732 
km, which is the distance corresponding to the experiment CERN-Gran Sasso[31], 
E = 50 GeV and in the case of LMA solution of the solar neutrino problems (LMA 
parameters: Am 2 2 = 10~ 4 eV 2 , tan#i2 = 1). For SMA, we used the parameters 
Am 2 2 = 10~ 6 eV 2 , tan# 12 = 7.5 x 10~ 4 [32, 33]. 

Xis ~ 0.056, 
x^ IA ~ 0.0019, 
x{™ A - 1.9 x 10~ 5 . 



1. v e — Vn channel 

In order to distinguish the characteristic signature of this kind of models, we first 
discuss the behaviour in different extreme cases of the CP asymmetries, Acp defined 
as follows 

P - P_ _ p - p— 

CP = p TpTT = P e , + P w (74) 

To get the analytical results, we should remember that we expand Pf^ 1 to second 



order in S13 and P^ p to first order in S13 and we assume that £21 "C £31 <C 1 =>- L <C 



(GeV) x 262 km. P^ is the probability of oscillation between v e — due to SM 
and P^ P is t ne probability of oscillation v e — due to New Physics. 

Let us consider the case when 613 = 0. In this case all CP violating phases are 
due to New Physics. As in ref. [10], we consider 2 different limits: one for large value 
of S13 (close to the experimental bound: \U e s\ < 0.16 [34]) and the other for small 
S13 (typically small means for instance, taking the LMA solution for solar neutrinos, 
s 13 <0.01). 

In the "large" S13 limit (221/^31 <C ^U^U^z) / (^2^2)!), the SM probability is 
given by 

p e T = ^ 2 31 \u e3 u; 3 \ 2 . (75) 

The New Physics CP asymmetry A^p then reads 5 



As in the limit <5 13 = 0, all elements of Umns are real, we simplify the notation omitting the * in 
the following. 
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/ e d* _|_ e s 

Aqp ~ — - Im ' Se e/J 



^31 \ ^3^3 



9 e4 9^ sm(5 u - 5 24 ) + 6> M 6> r4 sin5 2 4 (77^ 

#31 \\Ue3UfjzJ \U e3 
+^r4Sin5 14 ^jV (76) 

It is important to notice that A^p seems to be enhanced by factor. In fact, 
using eqs(68-69), it is easy to check that the U e % in denominator is cancelled once we 
replace the Umns elements in terms of the 9y angles. 

In the " small " s 13 limit (x 2 i/x 31 > \{U e3 U^) / (U&U^)]), the SM probability of 
oscillation is given by 

P c f ! = 4x 2 21 \U e2 U ll2 \ 2 . (77) 
For the CP asymmetry, we find 

1 f e d * + e s 

9 e4 9^ sin(5i4 - <5 24 ) + 9^9 tA sin5 2 4 



a?2i \\Ue2J \U, 



1 ni 



+e ^ Aah \€^))- (78) 

In case of bimaximal mixing, eqs. (76-78) can be computed making an expansion 
in S13 and using eqs(69): 

A^p ~ 1 a r4 s 13 sin</> + 0(sj 3 ). (79) 

2 p5M 



Both limits ("large" and "small" S13) are recovered using the appropriate P e/J 



SM 

From this equation, it is easy to get the limit of S13 = 0. 



A^is = 0) ~ « (80) 

with z = 2, 3. Thus we can conclude that CP asymmetries in the — v e channel are 
unobservable for s 13 = due to the constraint on the non-observation of rare muons 
flavour changing decays. 
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2. v T — Un channel 



Let us now consider CP violation in the v T — channel. In this case, the interesting 
limit is S13 -> A s c ^ = 0). The SM probability is 



Pl\ ~ 4sin^ 1 |^3| 2 |^3r, (81) 
leading to an induced asymmetry 

Acp * — (Y^K) 0e,9, A sin(5 14 - 5 24 ) + B^Ort sm6 24 ( %*) 

+9 ei 9 T4 smS u ( %± H . (82) 



We see that the last two terms are non-vanishing in the limit 812 ~ #24 and could 
then give observable effects at neutrino factories. Indeed, even for s 13 = 0, can 
be expressed using eqs(68) as 

A%£(s 13 = 0) ~ —a # r4 sin<5 2 4 (83) 

s 23 

C. Oscillations in matter 

In the general case, matter effects [35] have to be taken into account in neutrino 
oscillations. They are due to forward v — e scattering, the Lagrangian describing them 
coming from the following 4-fermion interactions 

C matt er = 2V2G F (5 ae 5 l3e + ) (V a ^P L e) {e-y„P L v p ) . (84) 

Using the notation given above, the different coefficients read 

(85) 
(86) 

4r£- (87) 

Thus we see that the first order corrections due to New Physics appear in e™, e™ and 
e™ while e^,e™ and e™ appear at higher order in New Physics interactions. 

■ 



m 
e e/3 


ANP _ d 

— ^efS — e ef3i 


m 


\/\NP 2 
| efi,er\ i 


m 


aNP aNP* 
— A ep, A eT — 



The general expression for the probability to have v s a — > vi is given by 



2 



P^ = |(^|^V^ 7 |z/ 7 )p. 
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Matter effects can be described using the following Hamiltonian 



H, 



a/3 



2E V 



(Umns) 




Amli 



{Umns% + b 







Am 



31 J 



jm jm ,m 



V 



eft 
in* 



with 6 = 1\p2G Ffi e E u and n e the electron density of the matter. 



IV. RESULTS 



/i/i fir 



[IT ^TT 



\ 1 



59) 



In this section, we shall numerically estimate the values of the CP asymmetries 
for the different cases studied in previous section, neglecting matter effects. We 
shall focus our analysis on the signal-to-noise ratio to emphasize on the possibility to 
measure such CP asymmetries in neutrino oscillations, in particular for the v e — 
channel, at future neutrino factories. For that, one proceeds as in ref. [36] defining the 
observable 

N[^}/N [e~}~N[ f i + }/N [e + } 



A 



where 

N\pT> + 
N [e-'+ 
with a 



CP 



N^NtE, 

7rL 2 m 2 
/* 

7rL 2 m 2 

H 1 



N[fi-}/N [e-}+N[fi+}/N [e- 



(90) 



12 



E u 



12 



E. 

E> 
E„ 



K 
K 

E„ 



<7cC,Cc( E ») P ^( E ») dE »> ( 91 ) 

cr C ncc(E u )dE u , (92) 



cc,cc\^ 



{E v 



CC,CC^v 



E v and a 



cc,cc 



given, respectively, by 0.67 x 



10~ 38 cm 2 /GeV, 0.34 x 10~ 38 cm 2 /GeV. iV M is the number of useful muon decays, N T 
is the number of protons in the target detector and E t h the threshold energy of the 
detector. 

The statistical error AAcp is given by 

1 



AA 



CP 



y/N\p.-] + N\n+] 



(93) 



The signal-to-noise ratio is given by Aqp/ AAcp- In order to illustrate the possibility 
to detect Acp at future neutrino factories, it is useful to compute the signal-to-noise 
ratio using the analytical results we have for Acp for the different cases studied in 
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previous sections. These results were obtained making an expansion in s 13 (expand- 
ing P^ 1 to second order in S13 and P^f to first order in S13) and were valid for 
short distances, £31 <C 1 =>• L <C E^(GeV)x262 km, (typically at neutrino factories, 
= 50GeV [4]). To get this number, we use Am^ = 3 x 1CT 3 eV 2 . Within these 
assumptions, A C p and P^ 1 have a very simple energy dependence 



A CP ~ £„, (94) 
^ M - ^- (95) 

Using this energy dependence, it is easy to integrate eqs. (91-92). The final signal-to- 
noise ratio reads 

To get an estimation, let us evaluate this expression for a total of 10 21 useful muons 
(Nfj) with an energy = 50 GeV and a 40kt detector (A^ « 1.1 x 10 33 ) and for 
= 732 km. 

^g- « 2 x itfAcPmjp™^) x (97) 

Usually, it is said that the signal could be distinguished from the background 
noise at 99% C.L. if the signal-to-noise ratio is bigger than three. Using eq.(79) and 
imposing the signal-to-noise ratio to be bigger than three, an lower bound on a and 
9 T 4 can be found, 

AcP > > 3 =>• a 6 t4 sin0 > —2 x 10~ 3 -^. (98) 



AAcpJ S13 Lgs 

As a 6> r 4 ~ Sf ~ /5 rT < 0.0013 due to lepton universality (see eqs(24)), the above 
lower bound is never satisfied for future neutrino factories as they are planned [4]. 

Thus, the constraints on FCNC in the leptonic sector and violation of lepton uni- 
versality make the New Physics contribution to CP-asymmetries in the v e — channel 
unobservable, in models with extra isosinglet charged leptons, with the present design 
of neutrino factories. 

We should stress at this point the importance of a measurement of CP violation 
in the — v T channel in neutrino oscillations. In this case CP violation mediated by 
the new charged lepton could be noticeable even in the limit S13 = what would be 
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a clear signature of physics beyond the SM. Indeed, proceeding in the same way than 
before, one gets 



\AA CP J ™ ^ '2 L 

This lower bound can be translated as a lower bound on S- 




(99) 



S 2 3 = Prr > \ x 10 



L 



(100) 



This lower bound is still marginally compatible with the upper bound on (3 TT coming 
from lepton universality. So, in principle, there is still a small window for observation 
of vector- like charged leptons effects on CP asymmetries in v^ — v T channel at neutrino 
factories. We should emphasise that as 5-factories produce as many tau pairs events 
as B — B events, the expected improvements on rare lepton flavour changing tau 
decays at BELLE experiment for instance should allow us to close this window in a 
very near future if no lepton flavour changing tau decays are observed. And inversely, 
the observation of a such events at a I?-factory will be a strong motivation to adapt 
the design of future neutrino factories to be able to probe CP asymmetries in — v T 
channel for the kind of model studied in this paper. 

V. CONCLUSION 

In this paper, we have studied models inspired by extra-dimensions where new 
particles naturally arise as isosinglet charged leptons. We have described a general 
perturbative approach to compute the effects of a vector-like charged lepton on neu- 
trino oscillations and their potential signature at future neutrino factories taking into 
account all constraints coming from FCNC and lepton universality. This approach is 
based on model-independent formalism introduced in ref. [10, 11] and can be applied 
to any extension of the SM with vector-like particles (neutrinos or charged- leptons) . 

We have shown that in case of leptonic maximal mixings, new physics effects in 
v e — CP asymmetries are significantly enhanced. But due to FCNC and lepton 
universality constraints, such effects are out of reach of current neutrino factory de- 
sign. In this class of models, we can expect to observe some signal from new physics 
at 1 — 3cr in — v T channel if rare flavour changing tau decays are seen in a very 
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near future at 5-factories as BELLE experiment. If such events are seen, there will 
be a strong motivation to slightly adapt the neutrino factory design to improve their 
sensitivity to CP asymmetries in — v T channel in order to test this class of models. 
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APPENDIX A: KALUZA-KLEIN EXPANSION OF BULK FERMIONS 

As we have emphasized in the Introduction, models with extra dimensions repre- 
sent a well motivated framework for physics beyond the SM with extra vector-like 
fermions. We review in this appendix the KK expansion of bulk fermions in the class 
of models introduced in [7]. Let us consider a five- dimensional model with the fifth 
dimension compactified on the orbifold S 1 /Z 2 , which is a circle of radius R with the 
Z2 identification y <-> — y or, equivalently, an interval < y < ttR with two bound- 
aries, the orbifold fixed points. A fermion in five dimensions is vector-like, the Dirac 
representation of 5"0(1,4) being irreducible, thus it admits a bare Dirac mass which, 
in order to be non-trivial, has to depend on the extra dimension. 

The integral in the fifth dimension of the five-dimensional Lagrangian results in 
the following four-dimensional Lagrangian (we use the "mostly minus" convention for 
the metric and 7 4 = i7 5 ) 

rR 



£ = I dy^ 

o-KR 



i^d N - M{y) 



dy 



f * R (d y - M(y)) V L + *l (-d y - M(y)) ^ R , (Al) 



where we have split the vector-like fermion into its two chiral components \P = ^ + 
defined by 7 5 ^l,_r = ' l,r- Note that in order to have a dynamical field in 
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the extra dimension, that is the term ^Ldy^R + h.c. 7^ 0, the two chiralities of a 
fermion necessarily have opposite Z2 parities to cancel the change of sign of d y . Thus 
without loss of generality we can choose the mass term to be odd M(—y) = —M(y). 
In particular we take it to have a multi-kink structure 

f M, < y < na, 
M(y) = { ~ y - (A2) 

I —M, na<y < nR, 

with < a < R. Hermiticity of the Lagrangian requires M to be real but it can 
be of either sign. To complete a four-dimensional description we expand the five- 
dimensional fields in a real complete orthonormal basis in the fifth dimension with 

the coefficients of the expansion being four-dimensional fields, the KK modes, 

00 

9 L>R (X, V) = -J= E fn '^lU^ ( A3 ) 

V 7rR n 

n=0 

Inserting this KK expansion in the Lagrangian we obtain the action of an infinite 
tower of four-dimensional vector-like fermions, except for the zero mode which will 
be chiral, 

C = J2 S^L + ^1 ~ m n [^^R + *? } *L ] ] } , (A4) 

n 

provided the following orthonormality and eigenvalue conditions are satisfied 

ttR cL eh rirR t-R i-R 

dy ^ = / dy ^ = 5 nim , (A5) 



TlR J txR 
[±d y -M(y)]ft R = -m n fK> L . (A6) 

The KK spectrum for this problem has been computed in [7] and can be summarised 
as follows. (In the following the upper (lower) sign stands for LH (RH) fields.) 

There is a massless zero mode for the even chirality with exponential localization 

fo' R (y) = A LtR exp[ T M\y - ira\], (A7) 
where the normalization constant is given by: 



.4 



±2MttR 



L,R 



2 - exp[=F2Mvra] - exp[=F2M7r(i? - a)] ' 

Notice that a LH zero mode is localized at the intermediate brane for M > and 
simultaneously localized at both orbifold fixed points with exponential suppression 
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in the intermediate brane (thus the designation we use hereafter "multi-localized") 
for M < 0, the opposite happens for a RH zero mode. The rest of the spectrum is 
vector-like with the first massive mode having distinct properties in the case that 



T2M7ra(R - a) > R. (A8) 

The condition (A8) coincides with the multi-localization of the zero mode (indeed we 
use this condition as a quantitative definition of multi-localization in this particular 
problem). In fact the multi-localization of the zero mode is intimately related to the 
special properties of the first massive mode. When the condition (A8) is satisfied, the 
even component of the first KK mode is also exponentially multi-localized 

ft R[cvcn \y) = { \ a A « A (A9) 



B (e-fcy + ggfe/Mv-^*)) , na < y < nR, 



where A and B are related by continuity of f( cvcn ^ at y = na provided the wave 
function does not vanish at this point and by continuity of the derivatives if the 
wave function vanishes at the intermediate brane. The parameter f3\ is the positive 
solution of the eigenvalue equation (with the upper sign valid for LH components and 
the lower one for RH fields) 



P 1 _ e -^« = T M 1 + e~ 2 ^ H - e" 2 ^ a - , (A10) 

for even fields and the corresponding with the change M — > —M for odd fields so 
that the two chiralities of a fermion KK mode have the same mass as they should. 
The mass of the first KK mode is mf = M 2 — /3f and is, provided condition (A8) is 
satisfied, always positive and smaller than M 2 . In particular it can be seen that in the 
case of strong multi-localization the mass of the first KK modes goes exponentially 
to zero 



2M 2 ^±2,Ma + e ±2«M(R-a) _ 2e ±2.A/fl] ^ Q frtf ^ ( AU ; 



It is therefore effectively decoupled from the compactification scale (thus from the rest 
of the KK spectrum) in this limit. In the case that the multi- localization condition is 
not satisfied this first state has the same properties as the rest of the spectrum which 
consists on oscillating (thus not localized) states with masses greater than M. 
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Let us finish this short review of the KK description of extra-dimensional theo- 
ries writing the Yukawa couplings. The Lagrangian involving the Yukawa couplings, 
taking the Higgs to live at the y = boundary, reads 

-C Yuk = / dy %){A (5) #x$ + h.c.} 
Jo 

\(5) 

= ^E $ {/» L (°)/^(°) + h - c -} . ( A12 ) 

nm 

where we have chosen \1/ L and \r t° be even fields and in the second equality we have 
expanded in KK modes. 
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